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Abstract 



'nI" ' It is well-known that f{R) theories in Einstein frame is conformally equivalent to 

■ quintessence models in which the scalar field minimally couples with gravity. If there 

exists a matter system in Jordan frame, then it interacts with the scalar field in Einstein 
qh| frame due to the conformal transformations. This interaction, in general, may lead to 

changes of fundamental constants. Here we will consider possible time variation of fine 
bX), structure constant in a general f{R) theory. We will use observational bounds on these 

variations and argue that it provides a criterion for constraining f{R) models. 

> , 

^ : 1 Introduction 
in 

! Recent observations on expansion history of the universe indicate that the universe is expe- 
2 I riencing a phase of accelerated expansion [1]. This can be interpreted as evidence either for 
Q I existence of some exotic matter components or for modification of the gravitational theory. In 
the first route of interpretation one can take a perfect fluid with a sufficiently negative pressure, 
dubbed dark energy [2], to produce the observed acceleration. In the second route, however, 
^ '. one attributes the accelerating expansion to a modification of general relativity. A particular 
^ I class of models that has recently drawn a significant amount of attention is the so-called f{R) 
gravity models [3] [4] . These models propose a modification of Einstein-Hilbert action so that 
the scalar curvature is replaced by some arbitrary function f{R)- 

It is well known that conformal transformations can be used to recast f{R) theories of gravity 
into the form of Einstein gravity together with a minimally coupled scalar field. The original 
variable is called Jordan conformal frame while the transformed set, whose dynamics is de- 
scribed by Einstein equations, is called Einstein conformal frame. If one formulates a f{R) 
theory with some matter systems in the Jordan frame, then the conformal transformation 
induces a coupling of the scalar field with the matter sector in the Einstein frame. One of the 
important implications of such a coupling is the possible changes of constants of nature. If the 
effective mass of the scalar field is sufficiently small the interaction of the field with ordinary 
matter results in the time variation of constants of nature over cosmological time scales [5]. 
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Among all possibilities, we restrict ourselves to time variation of fine structure constant as a 
consequence of interaction of the scalar field with electrodynamics. We intend here to use the 
current bounds on a running fine structure constant to constrain a general f{R) theory in the 
Einstein frame. 



2 Quintessence and time variation of fine structure con- 
stant 

Many years after Dirac's proposal [6] that fundamental constants may vary with time and/or 
space, there is now some theoretical frameworks which predict such changes. For instance, 
one of the most interesting low energy features of string theory is the possible presence of a 
massless scalar field, dilaton or moduli fields, whose vacuum expectation values define the size 
of the effective coupling constants. Independent of this framework, Bekenstein [7] has also 
formulated a theory for varying fine structure constant which has been recently generalized 
and apphed to a cosmological setting [8] [9] . Moreover, there is a large class of dark energy and 
quintessence models which invoke scalar fields with wavelengths comparable to the size of the 
universe. This scalar field may then interact with matter and radiation so that changes of its 
background value induce variation of coupling constants [5]. 

In recent years, there is an increasing interest to models concerning time variation of fine 
structure constant due to a non-vanishing coupling of quintessence to the electromagnetic 
field. In these models, one can combine cosmological data and terrestrial observations to place 
empirical constraints on the size of the coupling or even explore the conditions that the scalar 
field exhibit a tracking behavior [10] [11] [12]. The general form of the action in these models is 
the following ^ 

where R is the curvature scalar and is Lagrangian density of cosmic matter which we take 
it as a perfect fiuid with energy density Pm and equation of state parameter a;, and 

L<p = ld''<l>d,(l> - V{<j)) (2) 

L^p = -^A{<l^)F,,F^^ (3) 

A{(f)) is a dimensionless function-^ of and F^^, is the components of electromagnetic field 
tensor. For spatially fiat Friedmann-Robertson- Walker (FRW) spacetime, the governing field 
equations for the scale factor a{t) and the scalar field are 

H'^lip^ + l^' + Vm (4) 

^We work in the unit system in which h = c = SnG = 1 and our sign convention is (-+++). 
■I-We may call it a gauge function since evolution of fine structure constant strongly depends on the choice 
of this function. 
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H^-^{Ljp^ + (P^) (5) 

^ + 3i/0 + ^ = O (6) 

where H = -. These field equations depend on the choice of the functions V{(f)) and A{(f)). 
These functions characterize the evolution of the scalar field and the effective fine structure 
constant which is now given by 

with aQ being its present value and ^(^o) = 1- The choice of the gauge function A{(f)) is 
arbitrary although the choice of linear [11] and exponential [8] [12] functions are usual in the 
literature. It should give an appropriate temporal variations of fine structure constant which 
is consistent with observational constraints. We bring in the following two of these constraints 
which we will use in the next section^ : 

1) Estimates of the age of iron meteorites {z = 0.45) combined with a measurement of the 
ratio Re/Os resulting from the decay rate R^^^ Os^^^ gives [13] 

■ |< 10-^ (8) 



a 



2) We have also [14] 



a 



a 



< 4.2 X 10-^^ yr-^ (9) 



coming from comparing atomic clocks at present time = 0). The overdot denotes differenti- 
ation with respect to cosmic time. 

3 f{R) Gravity 

The action of f{R) gravity in Jordan frame in the presence of matter is given by 

S^lJ d^x^fiR) + J d'x^ L^{g^,, iP) (10) 

where f{R) is an arbitrary function of the scalar curvature. The Lagrangian density Lm 
corresponds to matter fields which are collectively denoted by ijj. We apply the following 
conformal transformation 

9,.. = ^9i.u, ^ = f'{R) (11) 

This together with a redefinition of the conformal factor in terms of a scalar field 

(h^^lnfl (12) 



§For a more complete list of these constraints see, e.g., [12] and references therein. 
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with P = ^Jl, yields [15] 

S=]^j d^x^g{R- r''V^<l>V.<j> - Vm + I d^'x^g e'^^'t' L^{g^,, ^) (13) 

In the action (13), is minimally coupled to gf^u and appear as a massive self-interacting scalar 
field with a potential 

V{<P) = le-'^Mm] - le-'^'^f{r[m]) (14) 

where the function r(f2) is the solution of the equation f'[r(n)]—n = [15]. Thus the variables 
{Qhui'P) provide the Einstein frame variables for f{R) theories. This Einstein frame represen- 
tation is equivalent to a quintessence model with a specific potential function which, due to 
(14), is closely related to the function f{R). 

It is important to note that in the action (13) the scalar field interacts with matter sector via 
the function e~^^'^ . If we take Lm to be the Lagrangian density of electromagnetic field, then 
the action (13) will become similar to (1) with a gauge function A{(j)) = e'^f^^ . Thus one can 
generally state that any j[K) theory in the Einstein frame is equivalent with a quintessence 
model in which the scalar field interacts with matter sector via an exponential gauge function. 
However, it should be pointed out that in the model (1) evolution of the scalar field depends 
on the choice of two arbitrary functions ^(0) and V^(0) while in the model (13) it depends 
only on V^(0) which in turn is characterized by the function f{R). It implies that depending 
on cosmological evolution of 0, or equivalently on functional form of f{R), the fine structure 
constant could have many possible histories during evolution of the universe. It is therefore 
worth studying what possibilities are allowed by the function f{R) within the available obser- 
vational constraints. 

To proceed furthur, we note that (7) is equivalent to a = ao e^^"^ . We would like to use (11) 
and (12) to write this relation in terms of scalar curvature rather than the scalar field, namely 
that a = ao f'^{R) where the prime denotes differentiation with respect to the argument. 
Thus the rate of change of a is given by 

a . f(R) , , 

In a spatially flat FRW spacetime R = 6{H + 2H'^) and therefore R = 6{H + AHH). In terms 
of deceleration parameter q = —ad/d^ = —(1 -|- H/H) and jerk j = a^d/d^ we will have 



It is also possible to obtain the relative change of a, 

Aa _ a{z) — ao 
a CKo 

which is equivalent to 



(17) 



r{R)u-i (18) 



a 
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In this relation a{z) and /'^(-R) \z indicate fine structure constant and /'^(-R) at redshift z. We 
can now use (16) and (18) to study the impact of the constraints (8) and (9) on a particular 
f{R) model. 

Let us first consider the model [16] 

f{R)^R + XR,A- (19) 

where A and n are positive parameters and Rc is of the order of the presently observed effective 
cosmo logical constant. Here and in the following we will take R^ = eH^ with e being a constant 
of order of unity and Hq is the Hubble constant. We substitute this f{R) function into (16) 
and (18) and obtain 

Arv 

— - nAx"-^(nAx'^-^ + 2) (21) 

a 

where x = R/Rc = -(1 — q)h^ with h = H[z)/Hq. Using the observational fact that Hq = 
72 ± 8 km^^s Mpc~^ = 7.4 x 10~^° yr~^ [17], then the experimental bounds (8) and (9) imply 



a n(n- l)Axr ^ 

- ^=0= 12(jo - go - 2) 
a 1 + uXxq 



o=12(jo-go-2) "'V; ^717-1 < 10-' (22) 

— |.=0 45= nAa;"-'(nAa;"-' + 2) < 10"^ (23) 

a 

where we have set £ = 1 and the subscript represents the present value of a quantity. It is 
now evident that the relations (22) and (23) act as algebraic constraints on the parameters A 
and n for given values of the quantities jo, qo, q{z — 0.45) and H{z — 0.45). 
To estimate these quantities, we use the gold data set for Jq and go which gives jo = 2.75lJ;fo 
and go = —0.86 ± 0.21 [18]. On the other hand, the Hubble parameter H and the deceleration 
parameter g are related by 

H{z) = Ho exp [ l\l + q{u))d\n{l + u)] (24) 

If a function q{z) is given, then we can find the Hubble parameter in the redshift z. Here we 
use a two-parametric reconstruction function characterizing q{z) [19] [20], 



1 qiz + g2 



where fitting the mo del to the gold data set gives gi = 1.47li:|2 and g2 = -1.46 ± 0.43 [20]. 
Using this in (24) yields [20] 



v2 



H{z) = Ho{l + zrexp[^+^^^^] (26) 

At 2; = 0.45, (25) and (26) give g = 0.12 and H — 1.28ifo. As an illustration, the changes 
of ^ and ^ in terms of the parameters A and n are plotted in Fig.l. The figure shows that 
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2.x 10-^" 4.x 10-'" 6.x 10"'" 8.x 10"'" l.xlO" 



Figure 1: The changes of - and — in terms of n and A for the model (19). The vertical 
axes are scaled according to the observational bound. The figures show that for larger values 
of n the slope of the lines increase. This implies that for moderate values of n, X should be 
extremely small. 



the bounds (22) and (23) are satisfied for small values of the parameters. For instance, for 
A to be of the order of unity one should have n < 10~^. A moderate value of n requires an 
extremely small value for A which makes the model (19) to be hardly distinguishable from 
Einstein gravity. 

Now we consider the models presented by Starobinsky [21] 



f{R) = R-XRc{l-[l + {^fr} (27) 



and Hu-Sawicki [22] 



R 



f(R) = R-XR,{-}Jk^} (28) 



which can be written in a unified form [23] 

f{R)=R-XR^{l-[l + {^Tr^^] (29) 

This parametrization corresponds to (27) and (28) with n = 2, and /3 = 1, respectively. To 
apply the experimental bounds, we use (16) and (18) and obtain 

-U=12(jo-go-2)^^^^ ' ^ < (30) 



Aa An , ^ ^ -x-i-irAn 

1 2=0.45 



,™ x^-\x^ + l)-^-\—x''-\x^ + l)-7-'-2} < 10-' (31) 
a p p 

We have used the same data set for the quantities jo, go, <l{z = 0.45) and H{z = 0.45) and 
plotted (30) and (31) in Fig. 2. The bounds satisfy in the region of the parameter space which 
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(c) (d) 




Figure 2: The changes of ^ and ^ for the Starobinsky (panels a, b) and Hu-Sawiki (panels 
c, d) models. The vertical axes are scaled according to the observational bound. In contrary 
to the model (19), in these models the larger values of n correspond to decrease of the slope 
of the lines. 



corresponds to larger values of n and smaller values of A. For a given value of n, the acceptable 
range of A is [0, Ai] with Ai being the maximum value of A which satisfies the experimental 
bounds. If we take larger values of n the bounds (30) and (31) allow a larger Ai which means 
that the whole range of acceptable A is larger. The figures show that for a given n{n > 1), Ai 
for the model (28) is smaller than that for the model (27). For n near unity, the acceptable 
values of A are very small for the both models. 



4 Conclusion 

A possible modification of gravity is replacing the curvature scalar in the Einstein-Hilbert 
action by a general f{R) function. It is well-known that these f{R) gravity models are confor- 
mally equivalent with a class of quintessence models in which gravity is minimally coupled with 
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a scalar field with some appropriate potential function. In this Einstein frame representation 
of f{R) models, conformal transformations induce a coupling of the scalar field with ordinary 
matter which may potentially lead to variation of constants of nature. 

Although both quintessence models and /(/?) theories in the Einstein frame can lead to changes 
of fundamental constants, there is a basic different between these two alternatives. In the 
quintessence models, the coupling of the scalar field with matter system is characterized by an 
arbitrary gauge function A{(f)). Cosmological evolution of is then determined by A{(f)) and 
some potential function of the scalar field V{(f)). On the other hand, we have shown that in 
f{R) theories the possible variation of constants only depends on the functional form of f{R). 
This is due to the fact that the gauge function is fixed by conformal transformations and is 
given by an exponential function and the potential of the scalar field is determined by f{R) 
itself. 

Based on these arguments we have introduced a criterion that can be used to probe the cosmo- 
logical viability of f{R) theories. In this criterion, we have considered a general f{R) theory 
in the Einstein frame and restrict ourselves to temporal variations of fine structure constant. 
We have shown that for any f{R) gravity model variations of fine structure constant can be 
written in terms of some cosmological parameters (such as g, j and H) in an appropriate 
redshift. In particular, the relations (16) and (18) allows us to apply the experimental bounds 
on the rate and the relative changes of a for any f{R) gravity model. We have then used the 
gold data set to estimate g, j and H at redshifts z = Q and z = 0.45 for applying to some 
particular f{R) functions. 

Specifically, we have applied the above procedure to models of Starobinsky and Hu-Sawicki. 
We have shown that these two-parameter models satisfy the experimental bounds in the re- 
gions of parameter space which correspond to larger values of n and smaller values of A. This 
result is more restrictive than the constraints coming from local gravity experiments which 
set a lower bound only on n [24] . In our analysis both the parameters n and A should satisfy 
the algebraic constraints (30) and (31) so that a particular value for one parameter affect the 
whole of range of validity of the other. Moreover, we have shown that for a given value of n, 
there is an upper bound on the allowed range of A (Ai). For n = 1, Ai will be of the order 
of 10~^ in both models. It is important to note that the local gravity analysis gives n > 0.5 
as the allowed range of n [24]. Although the outcome of our analysis is consistent with this 
result, it however emphasizes that this parameter should take values sufficiently larger than 
unity in order that the models keep an appropriate departure from Einstein gravity. 
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